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ABSTRACT 



Furukawa et al. 



(120091 ) reported the existence of a large mass of molecular 
gas associated with the super star cluster Westerlund 2 and the surrounding HII 
region RCW49, based on a strong morphological correspondence between NAN- 
TEN2 ^^C0(J=2-1) emission and Spitzer IRAC images of the HII region. We 
here present temperature and density distributions in the associated molecular 
gas at ~3.5 pc resolution, as derived from an LVG analysis of the ^^C0(J=2-1), 
^^CO(J=1-0) and ^'^C0(J=2-1) transitions. The kinetic temperature is as high 
as ~60-150 K within a projected distance of ~5-10 pc from Westerlund 2 and de- 
creases to as low as ~10 K away from the cluster. The high temperature provides 
robust verification that the molecular gas is indeed physically associated with the 
HII region, supporting Furukawa et al.'s conclusion. The derived temperature is 
also roughly consistent with theoretical calculations of photo dissociation regions 
(PDRs), while the low spatial resolution of the present study does not warrant 
a more detailed comparison with PDR models. We suggest that the molecular 
clouds presented here will serve as an ideal laboratory to test theories on PDRs 
in future higher resolution studies. 



Subject headings: ISM: clouds - open clusters and associations: individual (Westerlund 
2) - HII regions: individual (RCW 49) 
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1. Introduction 

The process of massive star formation is still poorly understood, partly because the 
majority of massive star forming regions are in highly confused and crowded environments. 
Study of the ISM around young stellar clusters can help to clarify the astrophysical activity 
in such regions, and can provide invaluable information on the earliest stages of massive 
star formation. 

Super star clusters ( SSCs) are mas sive clusters with stellar densities exceeding 10^ 



stars pc ^ in their cores (jjohnson 



20051 ). SSCs strongly impact the Galactic ecosystem by 
supplying the ISM with energy and momentum via strong UV radiation, stellar winds and 
supernova explosions. SSCs are very rare in the Galaxy and there are only several known 



to date, although other galaxies harbor more (e.g. 



Kornei fc McCrady 



2009. and references 



within). Three are locat ed in the Ga . 



and the Central cluster; 



Figer et al. 



actic center (the Arches cluster, the Quintuplet cluster 



1999 



Genzel et al. 



20031 ) and the other two in the 



Galactic disk (Westerlund 1 and Westerlund 2; hereafter Wd2). 

Wd2 is located at (/, fe) = (284.2, -0.33), close to the tangent of the Carina Arm. Of the 
five known SSCs in the Galaxy, it is the only one for which associated molecular clouds have 
been identified (Furukawa et al. 2009, see also Dame et al. 2007). The estimated age of the 



4x10^ to ~3xlO^M„ f Ascenso et al. 



2007 



Mq in the form of stars of 1<Mq < 120 Mq (lAscenso et al. 



Whitnev et al. 


20oi) 


(Ascenso et al. 


2OO7I; 



2004), with several thousand 



Rauw et al 



20071). This 



rich stellar population includes 12 O stars and 2 WR stars. 

Wd2 ionizes an extensive and luminous HII region, RCW49. This has recently been 



observed with the Infrare d Array Camera (IRA C 



of the GLIMPSE survey (IBenjamin et al. 



20031 ). 



on the Spitzer Space T elescope part 



Churchwell et al. 



(120041 ) present IRAC 



images of a dust emission and discuss polycyclic aromatic hydrocarbon (PAH) in the region. 



-4- 



which is likely excited by the strong ultraviolet /optical radiation from the central cluster 
Studies of IR point sources suggest tha t there are about 300 Young Stellar Objects (YSOs) 



more massive than 2.5 Mq in RCW49 (IWhitney et al 



2004h 



In addition, there are observations in high energy b ands in the region . RCW49 has 



been detected in X-ray s in the 0.2 to 4.5 keV range by 



Goldwurm et al. 



(119871 ). in the 0.1 to 2.4 keV range by 



Hertz fc Grindlavl ( 



Bellorii fc Mereghetti 



468 X-ray point sources have been detected with Chandra by 



extended X ray feature of mainly therma. 



Tsujimoto et al. 



984) and 



1994h. and 



2003). An 



Suzaku X ray satellite (IFujita et al. 



origin has also been detected toward Wd2 by the 
20091 ). Observations with the HESS TeV gamma-ray 
telescope reveal an ext ended TeV gamma-ray source, HESS J1023-575, toward RCW49 
2OO7I ). Emission scenarios include the inter action of cosmic r ay particles 



( lAharonian et al 



accelerated by stellar wind or supernova with molecular gas ( IFukui et al 



20091) 



Molecular ob servati o ns pr ovide an invaluable addition to this rich set of multi- 



wavelength data. 



Damel (120071 ) used ^2CO(J=1-0) survey data taken with the CfA 1.2m 



telescope at a resolution of 8.8' to identify an extended giant molecular cloud (GMC) at a 
velocity of ~11 km s~^ as a possible parent GMC of Wd2. More recently, Furukawa et al. 
(2009; hereafter Paper I) carried out ^^C0(J=2-1) observations at a resolution of 90" with 
the NANTEN2 4m telescope of Nagoya University. These authors present morphological 
evidence that the HII region is associated with two complexes of molecular gas, one in the 
velocity range —11 to 6 km s^^ and the other between 11 and 21 km s^^. These authors 
argued for the robust association of these clouds with Wd2, and concluded that a loose 
association with the more extended 11 km s~^ cloud of Dame (2007), while possible, is not 
well supported by the data. The kinematic dist ance to these molecul ar clouds is estimated 



3), 



to be 5.4lj^'4 kpc by using the rotation curve of iBrand fc Blitzl (119931 ). and the masses are 
estimated to be around 10^ Mq each. Based on the mass and spatio-velocity structure 
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of the gas, Paper I suggested that a colhsion between the two clouds may have triggered 
formation of the stellar cluster Wd2. 

In this paper we present an analysis of molecular data of the ^^C0(J=2-1), ^^CO(J=1-0) 
and ^^C0(J=2-1) emission lines. Section 2 details the observations themselves. Section 3 
presents maps, derives line ratios and carries out a Large Velocity Gradient (LVG) analysis 
in order to constrain temperatures and densities in the clouds. We find strong evidence 
from the derived temperature distributions to support the conclusions of Paper I. In section 
4 we compare our molecular gas temperatures with the results of simple PDR models, 
before summarizing our results in section 5. 



2. Observations 

This study makes use of new 13C0(J=2-1) observations, the recent ^2C0(J=2-1) data 
from Paper I and archiv al ^^CO(J=1-0) data from the NANTEN CO Galactic Plane Survey 



( jMizuno fc Fukui 



20041 ). Both J=2-l transitions were observed with the NANTEN 2 4m 
sub- mm/mm telescope of Nagoya University, which is installed at Pampa La Bola (4800 m 
above sea level), in northern Chile. 

Observations of the ^^C0(J=2-1) and ^"^00(1=2-1) lines were conducted in February 
2008 and December 2008, respectively. The backend was a 4 K cooled Nb SIS mixer 
receiver, and the single-side-band (SSB) system temperature was ~200 K in ^^C0(J=2-1) 
and ~100 K in ^^C0(J=2-1), including the atmosphere toward the zenith. We used 
two 2048-channel acousto-optical spectrometers (AOS), resulting in a bandwidth of 390 
km s~^ and a channel resolution of 0.19 km s^^ at 230 GHz. The telescope half-power 
beam width was 90" at 230 GHz. The pointing was checked regularly on Ori KL, 
(a,5)Bi95o = {5^ 32™ 47.0'', —5° 24' 22"), and the applied corrections were always smaller 
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than 20"; usually less than 10". The target region was observed between elevations of 
30 degrees and 60 degrees. All observations were conducted using an on-the-fly (OTF) 
mapping technique, and each spatial region was mapped several times in different scanning 
directions to reduce scanning effects. The final pixel size of the gridded data is 30". The 
effective integration time per pixel is ~2 s in ^^C0(J=2-1) and ~4 s in ^^C0(J=2-1), 
resulting in rms noise levels per 0.19 km s"^ channel of ~1.5 K and ~0.5 K, respectively. 
The standard source p Oph East, (0,5)^1950 = (16^ 29™ 20.9", -24° 22' 30"), was observed 
for intensity calibration. Calibration uncertainties of 12% in ^^C0(J=2-1) and 8% in 
^^C0(J=2-1) are estimated from the daily variation in its observed brightness temperature. 

The ^^CO(J=1-0) line was observed by position switching as part of the NANTEN 
Galactic Plane Survey. The beam size at 115 GHz was 2.6', and the data are gridded to a 
2' grid with a velocity resolution of 0.1 km s^^. The rms noise per channel is ~0.9 K, and 
calibration uncertainty is estimated at ~ 10%. 

3. Results 

3.1. Distributions 

Figure [1] shows the integrated intensity distributions of the three CO transitions, 
overlaid on a three-color IR image from the Spitzer/IRAC GLIMPSE survey. Following 
Paper I, we show the molecular gas distribution in velocity ranges of —11 to km s~^, 1 to 
9 km s~^ and 11 to 21 km s~^ and hereafter refer to emission in these ranges as the —4 
km s~^, 4 km s^^ and 16 km s~^ clouds, respectively. However, it should be noted that the 
4 km s^^ and —4 km s^^ clouds do not form completely distinct entities in the 3D data 
cubes, and that spectra in the —11 to 9 km s~^ range are sometimes complex and blended. 

As described in Paper I, the Northeastern half of the 16 km s~^ cloud is located toward 
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the central regions of RCW49. Notably, the Southeastern edge of the cloud, between 
{l,b) = (284.33,-0.28) and (284.27,-0.40), shows a strong correlation with a bright 
filamentary ridge of Spitzer emission ~ 2' Southeast of Wd2. The remainder of the cloud 
extends 0.2 degrees to the Southwest, extending fully outside the area of influence of the HII 
region. The 4 km s~^ cloud shows an excellent morphological correlation with the large-scale 
shape of the HII region. Emission at this velocity range is entirely absent within a radius of 
3' of Wd2, and at larger distances from the cluster the molecular gas fans out to the North 
and South, matching well with the IR-bright gas. The —4 km s~^ cloud distribution also 
matches well with the IR nebula and in particular we note that a prominent CO clump 
centered on (284.22, —0.29) appears to be associated with a n umber of IR point sour ces, 



one of which has been identified as a YSO of spectral type B2 (jChurchwell et al. 



20041 ). 



New to the present study is the addition of i3CO(J=2-l) data from NANTEN2, as well 
as archival ^^CO(J=1-0) data from the older NANTEN Galactic Plane Survey. All three 
lines follow very similar distributions at the present spatial resolution, with ^'^C0(J=2-1) 
emission detected towards the brighter regions of the ^^CO clouds. This can be seen in 
Figure [H 



3.2. Line Intensity Ratios 

We convolve the i2CO(J=2-l) and i3C0(J=2-l) data with a Gaussian of FWHM 2.6' in 
order to smooth the cubes to the resolution of the ^^CO(J=1-0) beam. We then derive peak 
main beam temperatures, Tmb (K), and linewidths for the relevant emission components 
at each ^^CO(J=1-0) pixel position by fitting either single or double Gaussian functions 
to the line profiles. Two ratios, i?2-i/i-o = T^bl^^ CO {J =2-1)] /Tmb[^^ CO {J =1-0)] and 
-R12/13 = Tmbl^"^ CO {J =2-1)] /Tmb[^^ CO {J =2-1)], are then calculated. The uncertainties on 
these ratios are around 15%, as calculated from the propagation of the intensity calibration 
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errors, under the assumption that these errors are uncorrelated. 

Of these ratios the former is of most immediate interest, being a measure of the level 
of excitation of the molecular gas. We find that in the 4 km s~^ cloud -R2-1/1-0 ranges from 
0.42 to 1.34 while the 16 km s~^ cloud shows values of between 0.36 and 1.09. In the —4 km 
s~^ cloud, for which all three lines were significantly detected at only five 2' pixels, ratios of 
between 0.80 and 1.34 are observed. The spatial distribution of -R2-1/1-0 is shown in panels 
(a) and (b) of Figure [2l It is striking that molecular gas at smaller projected distances 
from Wd2 generally shows higher ratios. In both the —4 and 4 kms~^ clouds, molecular gas 
spatially coincident with the HII region shows markedly enhanced ratios, and in particular 
the southern part of the 4 km s~^ cloud shows the highest ratios toward the cluster in the 
north but lower ones in the south. This suggests that gas closer to the cluster is rn o re hig hly 



excited, and is quite consistent with the association suggested by 



Furukawa et al. 



(I2OO9I) 



Panels (c) and (d) of Figure [2] show the spatial distribution of -Ri2/i3- This ratio ranges 
from 4.1 to 13.0 across the entire molecular gas complex. Smaller values indicate higher 
optical depths and, as expected, are generally found towards CO emission peaks where 
column densities are larger. 



3.3. LVG analysis 



We may explore the issues of 



velocity gradient (LVG) analysis (iGoldreich fc Kwan 



temperature and density more robustly by using a large 



19741 ). LVG analysis is a method of 



calculating a simplified rad iative transfer code based on the approximate escape probability 



formalism of ICastorl ( 1l970l ). We assume a spherically symmetric cloud of constant density 
and temperature with radial velocity distribution proportional to the radial distance from 
the cloud center. We then utilize the line ratios measured above to solve the equations of 
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statistical equilibrium for the fractional population of the lowest 40 rotational levels of the 
ground vibrational level of CO. The velocity gradient, dv/dr, is taken to be 2.6 km s~^ 
pc~^, as estimated from the typical linewidth in a single 3.5 pc cell. In order to evaluate the 
dependence of temperature on the adopted parameters we tested three values of dv/dr: 1.2, 
2.6 and 5.0 km s^^ pc^^, while keeping the other parameters fixed. We find that Tk varies 
by about a factor of 1.6 for a factor of 2 change in dv/dr, indicating that the results do not 
depend strongly on dv/dr. The CO fractional abundance X(CO) is takei i to be 5 x 10~^ 



(ISakamoto 



2004h . 



19931 ). and the CO/ CO abundance ratio is taken to be 75 (iGiisten &: Philipp 



Figure [3] shows solutions for kinetic temperature, T^, and number density, tiq, at 
selected positions, labelled A to F in Figure O In Figure H] we plot the spatial distributions 
of Tk and no over the entire extent of the clouds, only excluding points for which no 
^^C0(J=2-1) emission was detected. (In the case of the -4 km s~^ cloud ^^C0(J=2-1) is 
only significantly detected at seven pixels and we therefore exlude it from the following 
discussion.) It may be seen that for both the 16 km s^^ and 4 km s~^ clouds the kinetic 
temperature is enhanced toward the Spitzer IR nebula and at smaller projected distances 
from the cluster. This is especially true in the case of the 4 km s~^ cloud, which shows 
temperatures as high as 70—150 K on the side facing Wd2. Tk in the 16 km s~^ cloud is 
generally lower, reaching peak values of 30—70 K toward the IR nebula. However, these 
values are still significantly enhanced with respect to the typical temperatures of ~10 K 
observed in molecular clouds without a local heat source. Moreover, the 16 km s~^ cloud 
shows a clear decrease in Tk with projected distance from Wd2, decreasing to to 10 K 
beyond the IR nebula at distances of 15—20 pc from the central cluster. The temperature 
in the Southern part of 4 km s~^ cloud also drops slightly to 30—50 K beyond the IR nebula 
at projected distances of 10—20 pc from the cluster, while the Northern part of 4 km s~^ 
cloud located toward the IR nebula within 10 pc of Wd2 emains high in temperature. This 



- 10 - 



is illustrated in Figure [5], which shows plots of Tk as a function of projected distance from 
Wd2. This relationship supports the association between the molecular clouds and the 
stellar cluster suggested in Paper I, while the lower in the 16 km s~^ cloud compared to 
the 4 km s~^ cloud may suggest that it is located further from the cluster. We also note 
that the 4 km s~^ cloud follows more closely the shape of the IR nebula as it extends in the 
North-South direction, whereas the 16 km s~^ cloud is elongated to the Southwest, with 
only its Northeastern most portions coincident with the Spitzer emission. The extent of 
the IR emission in each cloud with increasing distance from Wd2 is indicated in Figure 5 
and explains why the high temperature regions in the 4 km s~^ cloud extend to larger radii 
than in the 16 km s~^ cloud. 

In general it is sensible to be cautious about the accuracy of the quantity Tk when 
the energy levels employed are not as high as the temperatures derived. In the present 
case the J=l and 2 levels of CO are at 5.5 K and 16 K above ground level, respectively. 
Nevertheless, the temperature ranges indicated in Figure 3 demonstrate that the analysis 
is accurate enough to conclude that the kinetic temperatures are significantly higher than 
10 K, the typical temperature of molecular gas without an extra heat source other than the 
general interstellar radiation field and cosmic ray sea. 

Number densities are typically ~3000 cm~^ in both clouds, with the Southern part 
of of the 4 km s~^ cloud showing higher densities than the Northern part. In particular, 
the ^^C0(J=2-1) integrated intensity peak has an estimated of ~8000 cm~'^. It is also 
notable that the 16 km s~^ cloud shows a slight density depression close to Wd2, suggesting 
gas dispersal by the stellar cluster. 

In PDR regions, the abundance of ^^CO may change by a factor up to 2.0 relative to 
^^CO in regions oi Ay less than ~1 magnitude, where chemistry becomes complicated by 



partial ionization (jVisser et al 



20091 ). The present clouds have significantly higher Ay, of 
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up to 10 mag, and the fraction of the positions with Ay less than 1.5 mag is only ~ 15%, 
suggesting that the effects of ^^CO abundance variation may not be important. In order to 
test the effects of the ^^CO/ ^'^CO ratio we examine cases in which this ratio varies between 
40 and 150, instead of the 75 assumed in the original analysis. We find that the derived 
temperature is still significantly high; for example, = 70 K at (Position A) in the -4 
km s~^ cloud falls in the range 45-90 K, significantly higher than 10 K. Considering the 
large values of Ay noted above, this temperature range gives conservative limits. Thus we 
argue that the effect of the PDR on ^^CO abundances does not alter the key result that 
temperatures towards the HII region are high, and the association of the clouds remains 
valid. 



4. Discussion 

4.1. Comparison with PDR Models 

RCW49 is a classic example of a photodissociation region (PDR), in which the dense 
molecular ISM and dust are being heated by a strong FUV field. It is therefore instructive 
to make a simple comparison of the PDR surface temperatures and the molecular gas 
temperatures obtained above. 



We make use of the standard model of 



Kaufman et al. 



(119991 ). which parameterizes the 



PDR surface temperature of a constant density cloud in terms of the number density of 
hydrogen nuclei, no, and the FUV (6 eV < E < 13.6 eV) field at the cloud surface, Gq. We 
take no as ~6xl0^ cm~^, based on the LVG analysis above, and estimate Go from the Far 
Infra Red (FIR) luminosity of the nebula. Here we utilize high-resolution IRAS 60 fim and 
100 yum data from the Infrared Processing and Analysis Center (IPAC). The filter properties 
of IRAS allow us to combine these data sets to obtain the FIR intensity, Ifir, between 42.5 
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/iin and 122.5 fim, via the following relation (Helou et al. 1988; Nakagawa et al. 1998): 



Ifir = 3.25 X 10"^ X F(60/xm) + 1.26 x 10"^ x F(100/im) 



T. 



Where F(60 /im) and -F(100 fim) are the 60 /im and 100 fim fluxes in MJy sr~^. Assuming 
that the FUV energy absorbed by the grains is reradiated in the FIR, we then estimate the 
FUV flux Go impinging onto the cloud surface from: 



Go = 47r X Ifir/1.6 x 10"^ 
Where Go is in units of the Habing Field, 1.6x10"^ ergs cm~^ s~^. Following 



(2) 



Kramer et al. 



( I2OO8I ). t his formula assumes that ph otons with E < 6 eV contribute about half of the total 



heating (iTielens fc HoUenbach 



19851 ). but al so that the tota. 



re-radiated bolometric IR 



2OOII ) . These factors approximately 



luminosity is a factor of two larger than Ipm (jPale et al. 
cancel. 

At the hottest location in the 4 km s~^ molecular cloud (Z,6) = (284.33, —0.366), 
F(60 fim) and F(100 /im) are 2.3x10^ and 1.7x10^ MJy sr~^, respectively . Go is therefore 



estima ted to be 7.5x10'^ ergs cm ^ s ^ sr ^. Comparing with Figured] of 



Kaufman et al. 



(Il999l ). the PDR surface temperature is thus estimated to be several hundred K. 



This result, although slightly higher than the derived LVG temperatures in the 
molecular gas, is nonetheless highly consistent when we consider the fact that the observed 
CO emission will not originate from the active cloud surface. Instead it is likely to originate 
from deeper, and therefore presumably colder, regions of the cloud, which have not yet been 
penetrated by the dissociating radiation. 

The above comparison is necessarily crude. Future high-resolution studies with 
instruments such as ALMA may resolve molecular gas emission down to the size scales 
as the structures observed by GLIMPSE, and pave the way to more detailed model 
comparisons. 
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4.2. Embedded YSOs in the -4km s ^ clump 



Finally, we shall also briefly mention the prominent clump in the -4 km s~^ cloud 
located at (Z, 6) = (284.23, —0.36). This clump shows hints of higher values of -R2-1/1-0 
and in its periphery than in its center (see Figure 2a and 4a), suggesting that its 
outer regions are heated by the UV photons of Wd2. The clump itself also appears to 
be associated with several internal YSOs, with a maxim um spectral type of around B2 



(jChurchwell fc Glimpse Team 



2005 



Whitney et al. 



2004 ) . The t otal luminos ity of a single 



ZAMS B2 star may be estimated from the evolutionary tracks of 



IbenI (Il965l ) to be around 



2.1x10^^ erg s ^. For comparison, the cooling rate of the molecular clump is estimated to 



Goldsmith fc Langer 



be Ix lO^'^^'^ '' ergs s^^, where we follow the cooling rate equation of 
(119781 ) and assume a kinetic temperature of 20 K, a density of 4x10^ cm~^ and approximate 
the clump as a sphere of radius 2.5 pc. This flgure is of the order of 1% of the estimated 
B2 star luminosity. Considering that the molecular feature is not bright in the Spitzer 
PAH bands, we suggest that the interior of the clump is unlikely to be primarily heated 
by the UV radiation from the central cluster, and that the embedded YSOs are more than 
sufficient to heat the gas to it's present temperature. It is both desirable and important to 
obtain future detailed molecular observations at higher resolution to explore its detailed 
properties in order to better understand the star formation occuring within it. 



5. Summary 

We have carried out an analysis of three CO transitions - ^^CO(J=1-0), ^^C0(J=2-1) 
and ^^C0(J=2-1) - toward Westerlund 2, and derived temperature and density distributions 
for the two associated molecular cloud complexes identified in Paper I. The main conclusions 
are listed below. These conclusions reinforce Furukawa et al.'s suggestion that the molecular 
gas is physically associated with Wd2 and RCW49, and therefore provides support for 
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their idea that a cloud-cloud collision may have triggered the formation of this remarkable 
cluster. 

1. The ratio of peak main beam tempertures in the ^^C0(J=2-1) and ^^CO(J=1-0) 
lines, i?2-i/i-o, ranges from 0.42 to 1.34, and the ratio in the ^3CO(J=2-l) and ^2CO(J=1-0) 
lines, i?i2/i3, ranges from 4.1 to 11.8. In both the 4 km s~^ and 16 km s~^ clouds, regions 
toward the IR nebula and closer to the cluster show enhanced values of i?2-i/i-o, indicating 
a higher level of excitement closer to the stellar heating source. -R12/13 ranges from 4.1 to 
13.0 across the entire complex of associated gas. 

2. We have carried out a Large Velocity Gradient (LGV) analysis to estimate 
temperatures and densities in the clouds. Kinetic temperatures derived are as high as from 
~30 to ~150 K toward the IR nebula within 10—15 pc of the cluster both in the 4 km s~^ 
cloud and in the 16 km s~^ cloud. These high temperatures require a local heat source, and 
provide robust verification of the physical association between the clouds and the HII region 
and stellar cluster, as suggested in Paper I. Densities are typically around 3000 cm~^, with 
peak values of around 8000 cm~^. 

3. The 16 km s~^ cloud, which extends beyond the IR nebula as far as 20 pc from 
Wd2, shows a clear decrease in with projected distance from the cluster. This, in 
addition to the high temperature toward the IR nebula, supports the association between 
it and the molecular cloud. The Southern part of the 4 km s^^ cloud also shows a similar 
temperature decrease with distance, with lower temperatures observed outside the IR 
nebula at a projected distance of 15—20 pc from the cluster. Similarly, the Northern part 
of the 4 km s~^ cloud, which is located toward the IR nebula and lies within 10-14 pc of 
Wd2, also shows generally high temperatures. The fact that the 16 km s-1 cloud shows 
lower temperatures than the 4 km s~^ cloud suggests that it may be located at a somewhat 
larger distance from the cluster. 
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4. We have made use of simple PDR model calculations to estimate a PDR surface 
temperature of several hundred K. This result is considered to be consistent with the 
LVG-based temperature estimates from the CO emission, which should originate from 
deeper and colder regions of the cloud. 

5. The embedded YSOs at (/,6) = (284.23, -0.36) provide more than sufficient energy to 
heat the -4km s~^ molecular clump in which they are located. 



6. Acknowledgments 

NANTEN2 is an international collaboration of ten universities, Nagoya University, 
Osaka Prefecture University, University of Cologne, University of Bonn, Seoul National 
University, University of Chile, University of New South Wales, Macquarie University, 
University of Sydney and Zurich Technical University. The work is financially supported 
by a Grant-in-Aid for Scientific Research (KAKENHI, No. 15071203, No. 21253003, 
and No. 20244014) from MEXT (the Ministry of Education, Culture, Sports, Science 
and Technology of Japan) and JSPS (Japan Soxiety for the Promotion of Science) as 
well as JSPS core-to-core program (No. 17004). We also acknowledge the support of the 
Mitsubishi Foundation and the Sumitomo Foundation. This research was supported by 
the Grant-in-Aid for Nagoya University Global COE Program, "Quest for Fundamental 
Principles in the Universe: from Particles to the Solar System and the Cosmos" , from 
MEXT. Also, the work makes use of archive data acquired with Spitzer Space Telescope 
and IRAS data gained with Infrared Processing and Analysis Center (IPAC). Spitzer is 
controlled by the Jet Propulsion Laboratory, California Institute of Technology under a 
contract with NASA. 



- 16- 

REFERENCES 

Aharonian, F., et al. 2007, A&A, 467, 1075 

Ascenso, J., Alves, J., Beletsky, Y., & Lago, M. T. V. T. 2007, A&A, 466, 137 
Belloni, T., & Mereghetti, S. 1994, A&A, 286, 935 
Benjamin, R. A., et al. 2003, PASP, 115, 953 
Brand, J., & Blitz, L. 1993, A&A, 275, 67 
Castor, J. I. 1970, MNRAS, 149, 111 

Churchwell, E., & Glimpse Team. 2005, in Revista Mexicana de Astronomia y Astrofisica, 
vol. 27, Vol. 23, Revista Mexicana de Astronomia y Astrofisica Conference Series, 
ed. S. Torres-Peimbert & G. MacAlpine, 53-59 

Churchwell, E., et al. 2004, ApJS, 154, 322 

Dale, D. A., Helou, G., Contursi, A., Silbermann, N. A., & Kolhatkar, S. 2001, ApJ, 549, 
215 

Dame, T. M. 2007, ApJ, 665, L163 

Figer, D. F., Kim, S. S., Morris, M., Serabyn, E., Rich, R. M., & McLean, I. S. 1999, ApJ, 
525, 750 

Fujita, Y., Hayashida, K., Takahashi, H., & Takahara, F. 2009, ArXiv e-prints 
Fukui, Y., et al. 2009, PASJ, 61, L23 

Furukawa, N., Dawson, J. R., Ohama, A., Kawamura, A., Mizuno, N., Onishi, T., & Fukui, 
Y. 2009, ApJ, 696, L115 



-17- 



Genzel, R., et al. 2003, ApJ, 594, 812 

Goldreich, P., & Kwan, J. 1974, ApJ, 189, 441 

Goldsmith, P. F., & Langer, W. D. 1978, ApJ, 222, 881 

Goldwurm, A., Caraveo, P. A., & Bignami, G. F. 1987, ApJ, 322, 349 

Giisten, R., & Philipp, S. D. 2004, in The Dense Interstellar Medium in Galaxies, ed. 
S. Pfalzner, C. Kramer, C. Staubmeier, & A. Heithausen, 253 

Hertz, P., & Grindlay, J. E. 1984, ApJ, 278, 137 

Iben, I. J. 1965, ApJ, 141, 993 

Johnson, K. E. 2005, in lAU Symposium, Vol. 227, Massive Star Birth: A Crossroads of 
Astrophysics, ed. R. Cesaroni, M. Felh, E. Churchwell, & M. Walmsley, 413-422 

Kaufman, M. J., Wolfire, M. G., HoUenbach, D. J., & Luhman, M. L. 1999, ApJ, 527, 795 

Kornei, K. A., & McCrady, N. 2009, ApJ, 697, 1180 

Kramer, C., et al. 2008, A&A, 477, 547 

Mizuno, A., & Fukui, Y. 2004, in Astronomical Society of the Pacific Conference Series, 
Vol. 317, Milky Way Surveys: The Structure and Evolution of our Galaxy, ed. 
D. Clemens, R. Shah, & T. Brainerd, 59 

Rauw, G., Manfroid, J., Cosset, E., Naze, Y., Sana, H., De Becker, M., Foellmi, C, & 
Moffat, A. F. J. 2007, A&A, 463, 981 

Sakamoto, S. 1993, PhD thesis. University of Tokyo 

Tielens, A. G. G. M., & HoUenbach, D. 1985, ApJ, 291, 722 



Tsujimoto, M., et al. 2007, ApJ, 665, 719 

Visser, R., van Dishoeck, E. F., & Black, J. H. 2009, A&A, 503, 323 
Whitney, B. A., et al. 2004, ApJS, 154, 315 



This manuscript was prepared with the AAS macros v5.2. 



- 19 - 




ZB^.'- 20^.3 HhLJ ZB^.i JOA^ 1^.' If .7 IBA. I 3H J ^ 29^.^ ■ 



GalodiC InfiUfd* Cdw*^ 

Fig. 1. — CO integrated intensity contours overlaid on Spitzer three color image of RCW49. 
The left column shows ^^CO(J=1-0) data, the middle column ^^C0(J=2-1) and the right 
column ^^C0(J=2-1). The top, middle and bottom rows show emission in velcity ranges 
of 11.0 to 20.9 km s^^, 1.2 to 8.7 km s^^ and —11.0 to 0.3 km s^^ respectively. Contour 
levels are 10 + 10 K km s"^ for i2cOfJ=l-0). 15 + 10 K km s^^ for ^2COfJ=2-l) and 2 + 
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Fig. 2. — Maps of the intensity ratios -R2-1/1-0 (top) and -R12/13 (bottom) overlaid with 
^^C0(J=2-1) integrated intensity contours from Figure [H Panels (a) and (c) show the 4 km 
s~^ cloud and panels (b) and (d) the 16 km s~^ cloud. The -4 km s~^ component is shown 
in the inset blue dotted box. The red cross marks the position of Wd2. Letters refers to the 
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Fig. 3. — Contour plots of the LVG solutions at points A to F, as marked in Figure [2l Solid 
and dashed lines mark solutions derived from the ratios -R2-1/1-0 and -R12/13 respectively. 
Here the outer lines indicate 15% calibration uncertainties, and shaded grey areas are where 
the two solutions coincide. Panels A-F assume X{CO) of 5xl0~^ and dw/drof 2.6 km s~^ 
pc^^. Panels Ea, Eb and Ec show the solutions for point E under the assumption of values 
of X/(du/dr) of 1.0x10"^, 2.6x10"^ and 4.0x10"^ (km s"^)"^ pc, respectively. 
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Fig. 4. — LVG-derived temperature (panels (a) and (b)) and density (panels (c) and (d)) 
distributions overlaid with ^^C0(J=2-1) integrated intensity contours from Figure [H Panels 
(a) and (c) show the 4 km s~^ cloud and panels (b) and (d) the 16 km s~^ cloud. The -4 km 
s~^ component is shown in the inset dotted box. The red cross marks the position of Wd2. 
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Fig. 5. — Plot of relation between LVG-derived kinetic temperature and projected distance 
from Wd2. Panel (a): red points show the 16 km s~^ cloud. Panel (b): Black and blue 
points show the Northern part of the 4km s"^ cloud and the Southern part of the 4 km 
cloud, respectively. Pink bars show the extent of the IR nebula RCW49. 



